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We present a theoretical formulation of competition among electromagnetically induced trans- 
parency (EIT) and Raman processes. The latter become important when the medium can no longer 
be considered to be dilute. Unlike the standard formulation of EIT, we consider all fields applied 
and generated as interacting with both the transitions of the A scheme. We solve Maxwell equations 
for the net generated field using a fast-Fourier-transform technique and obtain predictions for the 
probe, control and Raman fields. We show how the intensity of the probe field is depleted at higher 
atomic number densities due to the build up of multiple Raman fields. 

PACS numbers: 42.50.Gy, 42.65.-k 



Multilevel atomic system interacting with several elec- 
tromagnetic fields can give rise to variety of phenomena 
that depend on the strength and detunings of the fields. 
Often, the various processes compete with each other, 
whereby some processes are suppressed or interference 
between processes renders the medium transparent to the 
applied fields [J]. Well known examples include compe- 
tition between third-harmonic generation and multipho- 
ton ionization YM , and four- wave mixing and two-photon 
absorption 0, Very recently, Harada et al. demon- 
strated experimentally that stimulated Raman scatter- 
ing can disrupt electromagnetically induced transparency 
(EIT), where the incident probe beam is depleted and 
new fields are generated via Raman processes Q. The 
disruption of EIT is important to understand because it 
may degrade the performance of EIT-based applications, 
such as optical memories and buffers, and magnetome- 
ters. In this paper we present a theoretical formulation 
that enables us to study the competing EIT and various 
orders of Raman processes to all orders in the applied 
and generated fields. 

The standard treatment of EIT |(| is based on the 
scheme of Fig. Q] where the atoms in the state \c) in- 
teract with a probe field of frequency ui. A control field 
of frequency lo c interacts on the unoccupied transition 
| a) <-» | b). The probe and the control fields are tuned 
such that 



W — LU C — LUfjc 



(1) 



This results in no absorption of the probe field provided 
the coherence pb c has no decay. This treatment assumes 
that the frequency separation Wf, c is so large that the 
interaction of the control field uj c (probe field lu) with 
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FIG. 1: (Color online) A schematic diagram of a three-level 
atomic system with energy spacing huibc between two ground 
states |c) and |6). The control field with frequency ui c and 
probe field with frequency u act on the atomic transitions 
\a) <-> \b) and \a) <-> |c), respectively. 



the transition \a) <-> \b) (|o) <-> |c)) can be ignored. At 
high atomic number densities or for strong fields, this 
approximation no longer holds, which is the situation we 
consider here. 

At higher densities, Raman processes start becoming 
important [3, la] , such as those shown in the Fig for 
example. The Raman generation of the fields at oj c — cube, 
uj+LUbc can further lead to newer frequencies like lo c — 2ojb c - 
In order to account for the Raman processes, we write the 
electromagnetic field acting on both the transitions as 



E(t) = £{t)e 



(2) 



where £(t) denotes the net generated field. At the input 
face of the medium £{t) has two components to account 
for both control and probe fields 



(3) 



Under the Raman-resonance condition l|TJl. we expect 
£{t) to have the structure 



£(t)=^ {n) e 



inuJbct 



(4) 
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FIG. 4: (Color online) Amplitudes of different Fourier com- 
ponents of the net generated field as a function of the atomic 
density of the medium. The normalized propagation length 
a C/7 = 200 is equivalent to actual length of the medium 
L=7.13 cm with an atomic density of n — 10 10 atoms/cm 3 . 
The other parameters of the above graph are chosen as: in- 
put Rabi frequency fig = 0.57, £p/£c = 0.5, A c = 0.0, 
7 = 9.475 x 10 6 , T bc = 0.0, uj bc = IOO7 and A = 766.4 nm. 



FIG. 2: (Color online) Diagrammatic explanation of the (a) 
Stokes and (b) anti-Stokes processes. The intermediate state 
is denoted by \a). 




FIG. 3: (Color online) Three-level A system interacting with 
the space-time dependent field E(z,t) on both the optical 
transitions. 



Thus, £(~V gives the strength of the Stokes process of the 
Fig. Ufa); £( +2 ) gives the strength of the process of the 
Fig. H[b); and £^ +1 ' describes the changes in the probe 
field. For low atomic number densities, we expect the 
usual results and therefore £ w £ and £(°' « £ c . 

To calculate the net generated field £{t) for arbitrary 
atomic number density, we have to solve the coupled 
Maxwell and density matrix equations. We consider now 
the situation as shown schematically in the Fig. [3J The 
applied field E(z, t) couples the excited state \a) to both 
ground states \b) and |c) . Here 27's represents rates 



of spontaneous emission. In a frame rotating with the 
frequency lo c the density matrix equation for the atomic 
system are given by 

Paa = (Pba + Pea) - i&>£ (Pab + Pac) - ^IPaa , 

Pbb = i^Spab - iftePba + ^IPaa , 

p ab = -[27 -iA c ]p ab +iQs(Pbb - Paa) + i&£pcb , (5) 
Pac = -[27 - l(A c - UJ b c)]Pac + i^-SPbc + i^s(Pcc ~ Paa) , 
Pbc = -(Tbc + iu b c)Pbc + i^sPac - i^SPba , 

where the detuning A c and the space and time depen- 
dent Rabi frequency fig of the generated fields are defined 
by 

d £ 

A c ^uj c -u ac ; Qe(z,t) = —-. (6) 

For simplicity, we have assumed d ab = d ac = d. The 
elements p ac and p ab in the original frame can be obtained 
by multiplying the solution of Eqs. JSJ by e~ luJct . The 
induced polarization V is given by 

V=(dp ab + dp ac )e-^ t . (7) 

The Maxwell equations in the slowly varying envelope 
approximation lead to the following equation for the gen- 
erated field 

{-o7 + ^i)= l 2^ +p -^ (8) 




FIG. 5: (Color online) The spectral amplitudes of different 
fields are plotted against the atomic density of the medium. 
The normalized propagation length a(/j = 2000 is equivalent 
to actual length of the medium L=7.1 cm when the atomic 
density n = 10 11 atom/cm 3 . The other parameters are chosen 
as: input Rabi frequency Qe ~ O.57, £ p /£ c = 0.5, A c = 0.0, 
r 6c = 0.0 and u) bc = IOO7. 

where a is given by 

a = 3A 2 n 7 /47r, (9) 

and n is the atomic density. The coupled equations (JSJ 
and JSJ are solved in the moving coordinate system 

r = t--;C = z. (10) 
c 

We have numerically solved the coupled set of equa- 
tions when all the atoms are initially in the state |c) and 
when the fields at the input face of the medium are given 
by ©. We calculate £(/,t) and do a fast Fourier trans- 
form to obtain the different Fourier components of the 
field at the output face of the medium. This procedure 
enables us to find how the probe and control fields evolve 
and determine when the Raman processes become im- 
portant. In the simulations, we have used parameters 
that are appropriate for 39 K vapor Q to avoid numerical 
problems. In this situation, the spontaneous decay rate 
of the excited state \a) 47 — 3.79 x 10 7 rad/s and the 
wavelength for the ground state |c) to excited state \a) 
transition A = 766.4 nm. 

EIT Vs Raman Processes In this section, we present 
the results of numerical calculations. In Fig. 21 we show 
result for the low-density regime. In this region, we notice 
almost no change in the probe field and thus EIT dom- 
inates. It is also seen that the Raman processes slowly 
start to build up, leading to the drop in the control field 
amplitude. 



FIG. 6: (Color online) The amplitudes of hyper Raman com- 
ponents at frequencies (u c — 2ui, c ), (w + 2u>bc) as a function 
of the atomic density. All parameters are same as in Fig. JSJ . 



We next consider the high-density regime, as shown 
in Fig. [3] This is the region when multiple Raman pro- 
cesses build up significantly 9J. Our numerical results 
are in broad agreement with the observations of Harada 
et aZ.pJ where they observe the depletion of the probe 
field and the generation of the Stokes field at (uj c — u>b c ). 
In particular, we see in fig.[S]that the generation of radi- 
ation at (u> c — u>b c ) is very important and the probe beam 
is depleted. We also notice a new feature - the probe ex- 
hibits some oscillatory character before dying out. This 
oscillation is due to the fact that any population that 
is transferred to the state \b) can produce a field at the 
probe frequency via the Raman process. When this hap- 
pens, the control field amplitude falls. 

In Fig. [ft] we show the build up of several hyper-Raman 
processes. The effect of a buffer gas on the generated field 
is shown in Fig.0 where it is seen that the amplitude of 
the probe field depletes faster in the presence of buffer a 
gas. On comparison of Fig. [S]and Fig. d we see that the 
amplitudes of the probe field and the generated Raman 
field become equal at aC/7 = 272 (without buffer gas) 
and 84 (with buffer gas). This is in agreement with the 
observation in Ref . Q . 

In conclusion, we have investigated competition be- 
tween electromagnetically induced transparency and Ra- 
man processes in a A system due to the cross talk among 
the optical transitions. We have demonstrated that the 
EIT-induced probe spectrum is very pronounced in com- 
parison to the higher order Raman sidebands for a low 
atomic number density. However, the generated Raman 
fields become dominant for an atomic number density 
that is only ten times higher. 
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FIG. 7: (Color online) The amplitudes of Stokes, control, 
probe and anti- Stokes fields as a function of the atomic num- 
ber density in the presence of a buffer gas. Here, we have 
scaled the amplitudes of probe, Stokes and anti-Stokes fields 
by a factor of two. The other parameters are chosen as: den- 
sity n = 10 10 atom/cm 3 , input Rabi frequency Q,£ — 0.57, 
Bp/Ec = 0.5, A c = 0.0, T 6c = O.OI7 and w bc = IOO7. A 
nonzero value of Tb c accounts for the buffer gas. 
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